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ABSTRACT 

Our aim is to study the evolution of the orbit of a star under the influence 
of interactions with an accretion disc in an AGN. The model considered 
consists of a low-mass compact object orbiting a supermassive black hole 
and colliding periodically with the accretion disc. Approximate calculations 
based mostly on Newtonian theory of gravity have been carried out by several 
authors to estimate the effects of circularization of initially eccentric orbits 
and their dragging to the disc plane. Here, we present the first step to a more 
adequate general rclativistic approach in which the gravitational field of the 
nucleus is described by the Kerr metric. The star is assumed to move along 
a geodesic arc between successive interactions with an equatorial accretion 
disc. We solved relevant formulae for the geodesic motion in terms of elliptic 
integrals and constructed a fast numerical code which, after specifying details 
of the star-disc interaction, enables us to follow the trajectory of the star for 
many revolutions and study the evolution of its eccentricity and inclination 
with respect to the disc. Lense-Thirring precession of the orbit is potentially 
a very important effect for observational confirmation of the presence of a 
rotating black hole in the nucleus. Our approach takes effects of the Lense 
Thirring precession into account with no approximation. 

The model of star-disc interactions has been suggested to explain the X- 
ray variabihty observed in the Seyfert galaxy NGC 6814. We briefly discuss 
on this subject. 

Key words: general relativity — black holes — galaxies: active — galaxies: 
individual (NGC 6814) 

1 INTRODUCTION AND MOTIVATION 

General interest in studying star-disc interactions in the nuclei of galaxies 
has greatly increased in the last years. This is partly due to the fact that 
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they appear important in explaining the X-ray variabihty of active galactic 
nuclei (AGN). Although it is generally believed that many galaxies, and ac- 
tive galaxies in particular, harbour massive black holes in their cores, there 
is no direct observational confirmation for this paradigm. The origin of this 
difficulty is apparent: complicated plasma physics of the matter swirling 
around the black hole makes it difficult to distinguish the effects of general 
relativity — although they may be essential for the mechanism of energy gen- 
eration itself. Potentially very important observable is the X-ray data on 
variability of active galactic nuclei (for the recent review see Wallinder, Kato 
& Abramowicz 1992). Although our understanding of the origin of X-rays 
is not satisfactory, it is often accepted that they are generated in the inner 
regions of the accretion disc where general relativistic effects are significant. 
When trying to describe these effects it is useful to separate details of the 
mechanism generating radiation (described in the local frame co-moving with 
the matter), i.e. the local physics of the interaction on the one hand, and 
observable effects as seen by a distant observer on the other. In our previous 
work (Karas, Vokrouhlicky & Polnarev 1992, Paper I) we developed a code 
which can be used in many astrophysically relevant situations to calculate 
images of various effects occurring in the close vicinity of the rotating (Kerr) 
black hole. Our code deals efficiently with problems treated originally by 
Cunningham & Bardeen (1973). All relativistic effects on photons (like grav- 
itational and Doppler shift of frequency and bending of light rays) were taken 
into account. As an example, we apphed the code to the "hot spot" model 
of the Seyfert galaxy NGC 6814 which tries to explain its X-ray variability 
on the scale of approximately 3 hours in terms of a bright orbiting spot (or 
spots) located on the accretion disc (Abramowicz et al. 1992), and we also 
discussed the case of a large number of spots with different intrinsic charac- 
teristics which may be relevant in explaining the X-ray variability of AGN 
on still shorter time scales (Abramowicz et al. 1991). 

In this paper we present a method for calculating the evolution of the 
orbital parameters of a compact star orbiting a massive black hole. Such a 
star may come from a binary or a cluster tidally disrupted by the central black 
hole (Hills 1988, Novikov, Pethick & Polnarev 1992). It can be deposited in 
a tightly bound orbit with a close pericentre where relativistic effects are 
important. The star interacts with the accretion disc only at the moment 
when it crosses the equatorial plane and this interaction weakly affects its 
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motion. Cumulative effects of successive tiny interactions circularize the 
trajectory and change the orbital plane into the plane of the disc. Mutual 
star-disc interactions have been discussed by Syer, Clarke & Rees (1991) as a 
possible origin of fueling and variability of AGN. Particularly, they discussed 
relative time-scales for circularization of the orbit and its "grinding" to the 
disc plane, and the final radius of the embedded orbit. Naturally, quantitative 
estimates depend on poorly known details of the interaction (cf. Zurek, 
Siemiginowska & Colgate 1992). A compact star cohding with an accretion 
disc is one of viable models for NGC 6814 (Abramowicz 1992, Sikora & 
Begelman 1992, Rees 1993). All present models possess their advantages 
and difficulties when compared to observational data and they need to be 
investigated in greater detail. General relativistic precession of orbital nodes 
(Lense-Thirring effect) — if detected — would strongly support models of AGN 
with rotating supermassive black holes. Lense-Thirring precession affects the 
inclined trajectory of a star, but it has been discussed only in special cases 
of free orbits with a large radius compared to the gravitational radius of 
the central black hole (Lense & Thirring 1918) and spherical orbits around 
an extreme Kerr black hole (Wilkins 1972). We consider the more general 
case of eccentric orbits which interact with the equatorial disc and we do not 
assume any particular value of the angular momentum of the central black 
hole. 

Star-disc interactions are a complex problem. Because we assume that 
the disc produces only a weak perturbation of free motion of the star, we 
can attack the problem in several steps. The present paper concentrates 
on an effective method for calculating the free motion of the star between 
its successive interactions with the disc. Simple examples of how star-disc 
interactions may change the picture are given in the last Section and they will 
be discussed in a forthcoming paper in greater detail. Detailed calculation of 
the Lense-Thirring frequency relevant for this model will be given elsewhere 
(Karas & Vokrouhhcky 1993). Once the motion of the star and its interaction 
with the disc are specified, we can apply the method of Paper I to compute 
the shape of the light curve or resultant spectrum. [See also Cunningham & 
Bardeen (1973), de Fehce, Nobih & Calvani (1974), Luminet (1978) and 
Laor & Netzer (1989)]. Recently, Fabian et al. (1989), Kojima (1991) and 
Laor (1991) applied analogous approach to study line profiles from accretion 
discs. For a more complete hst of references cf. Paper I. 
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Let us briefly describe the configuration of the modeL We consider a 
low-mass compact object (white dwarf, neutron star or black hole) orbiting 
the central massive black hole of the AGN. In our approach, we restrict our- 
selves to the assumption that the orbiting object moves along a time-like 
geodesic in the unperturbed background Kerr metric outside the equatorial 
plane. Thus we implicitly assume that (i) the orbiting object is sufficiently 
compact and/ or far from the central black hole (we neglect any coupling of 
the star's higher multipoles to the background curvature), and (ii) its mass 
is very small compared to the mass of the central black hole (we do not con- 
sider perturbations of the background metric). We also neglect the infiuence 
of radiating gravitational waves. In several astrophysical situations such 
assumptions may not be appropriate — see, e.g.. Carter & Luminet (1983), 
Luminet & Marck (1985), and Carter (1992) who studied tidal squeezing of 
the stars by the nearby black hole. Hartle & Thorne (1985) and Suen (1986) 
developed a scheme for multipole-tidal interactions of relativistic objects. As 
for the motion of close and comparable in mass black holes, see, e.g., D'Eath 
(1975a,b). We do not consider such extreme situations in this paper. Kates 
(1980) has shown that the star will move closely along a geodesic in the un- 
perturbed background metric for sufficiently long time provided the ratio of 
its mass and the characteristic reference length of the background metric is a 
small parameter. We have in mind situations where this parameter is of the 
order 10~^ or even less. This is important to note because we will attempt 
to follow the trajectory of the orbiting object for long periods of time. Next, 
we assume that the accretion disc is in the equatorial plane of the central 
black hole. We exclude thick disc models from our present considerations. 
Each time the object crosses the equatorial plane it interacts with the disc 
(provided the intersection is between the outer and the inner edges of the 
disc). In other words, we assume that the trajectory consists of arcs of free 
geodesic motion above (or below) the disc plane and impulsive changes of the 
orbital parameters at the moment of passages through the disc. The whole 
"physics" of the problem is compressed into the prescription governing the 
changes of the orbital parameters when crossing the equatorial plane. The 
interaction is assumed very weak which implies that relative changes of en- 
ergy, angular momentum and other quantities characterizing the orbit of the 
star are much less than unity in each single event. We should note that the 
geodesic motion in the Kerr space-time is integrable; thus, we do not ex- 
pect strong dependence of the shape of the orbit on initial conditions which 
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is typical for a chaotic motion. An alternative approach which employs a 
statistical description with appropriately averaged quantities is under prepa- 
ration. The separation of dynamical and physical aspects, as we introduce it 
in the present paper, appears very advantageous and it allows us to employ 
a fast method for computing the evolution of orbital parameters. 

2 THE MAPPING - DETAILS OF THE CALCULA- 
TION 

We consider the geodesic motion of a test particle (representing a star or a 
low-mass black hole) in the fixed Kerr background metric. We are interested 

only in short arcs of geodesic motion with the following boundary conditions: 
the initial point (indexed "i") lies in the equatorial plane (^^ = 7r/2 in Boyer- 
Lindquist coordinates), the final point (indexed "/") is the nearest successive 
intersection of the orbit with the equatorial plane. We employ a 'mapping' 
by which we understand an analytical algorithm to evaluate the final position 
{rf,(j)f) from the initial position {ri,(j)i) with constants of motion assumed to 
be given. In applications, we also need to know the transformation from 
initial to final velocities dr/dt, d6/dt, d(j)/dt to obtain the full starting infor- 
mation for the physical model of the interaction of the orbiting object with 
the accretion disc. The whole procedure is trivial in principle because the 
geodesic motion is separable in the Kerr background metric and the equa- 
tions of motion can be reduced to a set of ordinary first order differential 
equations (Carter 1968). Our main task is to handle the problem efficiently. 

Carter's equations involve squares of the velocities. As the star crosses 
the equatorial plane, the latitudinal velocity changes its sign periodically. In 
order to treat the case of the radial velocity, we introduce the sign function 
77 = sgn{dr/dt). Thus, our mapping is the analytical transformation 

We will also be interested in analytical evaluation of the delay in coordinate 
time which is necessary to pass from the initial to the final configuration in 
the disc plane: {tf — ti). In particular, this is important for reconstruction of 
the AGN photometric curve, provided similar periodic process as described 
here (star-disc interactions) arise its variability (Karas & Vokrouhlicky 1993). 
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It is worth noting that the code based on this mapping technique is op- 
timized as concerns both the speed and accuracy. The effective step of the 
method is the whole orbital arc and it cannot be made greater in principle. 
Moreover, the exact analytical solution of the problem is chosen as a sample 
function covering one integration step (instead of, for example, polynomi- 
als in Runge-Kutta methods). We remark that the name 'mapping' comes 
from analogous methods developed in celestial mechanics (e.g. Wisdom 1982, 
Murray 1986). 

We use the standard notation for the Kerr metric (Bardeen 1973). Quan- 
tities with the dimension of length in geometrized units are divided by the 
mass of the central black hole M, and they are thus made dimensionless. 
Time-like geodesies in the Kerr space-time can be integrated in the form 
(Carter 1968) 



tf - ti 
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The constants of motion, £ — —pt, ^ = p^j), and /C, can be expressed in 
terms of components of the four-momentum in the locally non-rotating frame 
(LNRF): 
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/C=[($-a£)2 + E(/f] (7) 

Components of the four-momentum in terms of direction cosines in the local 
sky of an observer at rest with respect to LNRF are 

= 7, = IV cos a, 

= 7i'sinQ;cos/9, p''^ = 7i'sinQ;sin/3, (8) 
where v is the tetrad velocity of the particle in LNRF and the Lorentz factor 

_ 1 

Specification of the initial conditions as a result of the star-disc interaction in 
the local frame co-moving with the matter of the disc requires another boost 
to the disc co-rotating frame (DCF). In the case of a Keplerian thin disc in 
the equatorial plane the linear velocity of the DCF with respect to LNRF is 

- (r3/2 + a)A ■ 

It is known that integrals (l)-(3) can be reduced to standard elliptic integrals, 
but only the simplest cases have been discussed in the literature. The explicit 
form of relevant formulae depends on the value of the constants of motion 
and initial conditions. We shall restrict ourselves to the most interesting 
astrophysical case: stable, energetically bound trajectories which cross the 
equatorial plane many times repeatedly. Thus, we assume < £ < 1. We 
exclude singular cases of orbits lying exactly in the equatorial plane (Q = 0; 
Bardeen, Press & Teukolsky 1972) or those intersecting the rotation axis 
($ = 0; Stoghianidis & Tsoubelis 1987); we also exclude the case of the 
extremely rotating black hole, a = 1. In the Schwarzschild case (a = 0) the 
geodesic is always planar, while a 7^ leads to the Lense-Thirring precession 
of the orbit and we have to take into consideration the dragging of the nodes 
[Wilkins (1972) discussed this effect for spherical orbits in the extreme Kerr 
case] . 

First, we simplify the integrals (l)-(3) to a form which can be directly 
found in standard tables of elliptic integrals (Byrd & Friedman 1971, Grad- 
shteyn & Ryzhik 1980, Grobner & Hofreiter 1965). For this purpose we need 
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to find the roots of R{r) and 0(^^). Real roots of are the turning points of the 
radial and latitudinal motion, respectively. We specify the initial point of the 
geodesic as (ij, rj, 7r/2, 0j) and we look for the final position {tf,rf,n/2,(f)f) 
which is in the equatorial plane again. The polynomial R{r) governing the 
radial motion is of the fourth order with -R(rj) > 0. In our case, R{r) < for 
r oo and r = 0. Thus we can find two real roots, G (0, r^), G (rj, oo), 
and for the remaining two roots we obtain the quadratic equation 

(£2 _ ^ [(^2 _ + + 2l r - ^ = 0. (9) 

Supposing all roots are real we can denote them, in a descending sequence, as 
Ti > r2 > r3 > r4. We exclude the possibilities of multiple roots because such 
situations are singular in the sense that they occur for precisely arranged 
values of E and The probability that interaction with the accretion disc 
will lead to such values is zero (in the measure sense). Thus we have 



_ /■ dr L r dr . , 

( r T, I I T Tr. I I T I I T — r ,\ 



In the case of two real and two complex roots (ri > r2 and rs, ra, respec- 
tively), we obtain 

/.= / / "^r (11) 

^/l-£''J {r^-r){r-r2)^{r-Xi? + X2 

where 

Xi = 3f^e(r3), X2 = ^^(rg). 

The latitudinal motion is governed by the polynomial Q^ilJ') = sm^ 9 0(6*) 
[Eq. (5)]. Solving the bi-quadratic equation in = cos^, 

a^l - S^)i^^ - [Q + 0^(1 - + $2j ^2 ^ g ^ (^2) 

we obtain the roots //+>//_> 0; The latitudinal motion is only possible in 
the region /i G {—ii-,fj,-) and 

T = f ^ 1 [ d^i 
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Analogously, the azimuthal motion can be solved in the form 
(f)f-(f)i^ [2{a£ - ^)A+ + ^B+] I+ + [2{a£ - + /_ + $J^ , (14) 

where 

^* " / (r - r±)fl(r)l/2 ' •'>' = j (I- ^2)e„(,i,)i/2 ' 



r± = 1 ± Vl - a2 . 



Finally, for the time coordinate we obtain 

tj -U = £ {Jr + 2Kr) + 2 [B+r+£ + a{a£ - $)A+] /+ 

+2 [B_r_£ + a{a£ - $)A_] /_ + 4£Ir + a^£Kf, (16) 



with 

W^' ^' = Je^' J e^(/i)V2- (1^) 

It is straightforward but moderately tedious to derive the explicit form of 
the mapping. We give relevant formulae in Appendix. 



3 THE EVOLUTION OF ORBITS - SIMPLE EXAM- 
PLES 

There are several interesting issues closely related to our problem which have 
not been fully understood as yet. In this Section we present simple examples 
of secular changes of orbital parameters of a star orbiting a supermassive 
(10^-10^ M0) black hole and interacting with the accretion disc. In partic- 
ular, wc study the eccentricity and inclination of the orbit, as they are in- 
troduced in Appendix. Differences from previous estimates which have been 
made within the framework of Newtonian gravity are found to be significant 
for nearly unstable orbits. These differences can be very subtle and inter- 
fere with the details of star-disc interaction and thus at first we adopted an 
extremely simplified (and perhaps unrealistic) description of the interaction. 
The Lense-Thirring precession of the orbit was included fully relativistically 
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with no approximation. This is important for our particular example we 
discuss below, NGC 6814. Mittaz & Branduardi-Raymont (1989) and Done 
et al. (1993) determined the periodicity in modulation of the X-ray emis- 
sion from NGC 6814 to be ~ 12, 000 seconds. Somewhat uncertain value 
for the mass of the central black hole was estimated to ~ 10^ Mq and thus 
the radius of the corresponding orbit is a few tens of the gravitational radius 
of the black hole — rather close to the horizon where approximation methods 
for the Lense-Thirring precession are no longer satisfactory. We should note, 
for completeness, that if the star-disc interaction is switched off, the eccen- 
tricity of the orbit remains constant in time and the value of the precession 
is exactly that of the Lense-Thirring precession. In the following examples 
we tune the strength of the star-disc interaction in such a way that the rel- 
ative change of orbital parameters is of the order f« 10~^ in each interaction 
and we follow ~ 10^ revolutions. Each of the following Figures shows the 
sequence of radial coordinates r for successive intersections of the trajectory 
with the disc — one point corresponds to one intersection, two intersections 
correspond to one revolution of the orbiter. (Alternatively, instead of the 
number of intersections N we could use coordinate time t for labelling the x- 
axis; Figures remain very similar in shape but the second possibility appears 
more adequate for plotting computed light curves which should be related to 
the observer's time at infinity.) Upper and lower boundaries of the distribu- 
tion of intersections are current values of the apocentre and the pericentre, 
respectively. Intersections are (seemingly randomly) scattered in the whole 
range between these boundaries due to the shift of the pericentre and Lense- 
Thirring precession (if a ^ 0). Frequencies corresponding to both of these 
effects are quantitatively studied in Karas & Vokrouhlicky (1993). 

In the first example we assume that, as a result of the interaction, the 
difference in the azimuthal components of the star and the disc material 
(evaluated in DCF) is decreased: 

^V^^ DCF > a* 0, DCF , (18) 

where a* is a phenomenological parameter ^ 1. Initially retrograde orbits 
(those with I > 90°) decrease the absolute magnitude of their $ component 

of angular momentum due to interactions with the disc and they either get 
captured by the black hole or become prograde. Once they are prograde, 
the star acquires angular momentum from the disc and moves away from 
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the black hole. Simultaneously, both the eccentricity and the inclination 
decrease, and the orbital period increases. The effect of energy dissipation 
due to crashing through the disc is not considered in Eq. (18). The model is 
of course inadequate when the inclination reaches zero and the star becomes a 
part of the disc. Figure 1 is an example of such an orbit. As mentioned above, 
the unperturbed geodesic motion in the Kerr metric is integrablc and thus we 
do not expect any dependence of the characteristic time for circularization 
or changing of the trajectory to the disc plane on particular values of the 
initial position or directions of velocity of the star. 

Our second example is a modification of the model studied by Syer et 
al. (1991). It is complementary to the previous one because energy dissi- 
pation is now considered while the difference in azimuthal velocities of the 
star relative to the disc material is ignored. We supposed that the star hits 
the disc supersonically and pulls some amount of the material with a mass 
Am ^ pdisc^disc^effsin"^/ out of the disc; here, pdisc and /idisc are the lo- 
cal density and thickness of the disc, respectively, and A^s is the effective 
cross-section for the star-disc interaction. The energy dissipated during the 
interaction is proportional to the kinetic energy acquired by the disc mate- 
rial, A£^diss oc Am(7DCF — !)• We assumed that the acceleration of the star 
which results from this interaction is anti-parallel to the velocity of the star 
and the corresponding change of the velocity is 



(Wc should note that the last equation for the star-disc drag becomes inap- 
propriate and must be modified when the motion of the star is subsonic and 
the disc material is directly accreted onto the star. The motion is highly 
supersonic with the Mach number of the order 10^ — 10"^ under conditions we 
consider.) Naturally, pdisc and /idisc depend on the disc model. Because we 
do not want to enter into these additional details we assumed that they, as 
well as Aeff, are constant. [We have also carried out computations using the 
density profiles corresponding to the Novikov & Thorne (1973) relativistic 
thin disc model which yielded only moderate modifications to the results.] 
Figure 2 illustrates two typical cases — both orbits are initially prograde with 
(a) / = 35° and (b) / = 80°. In general, the final radius of the orbit can be 
either larger (for small values of the initial inclination) or smaller (for large 



Ai;= 



oc 



Pdisc ^disc ^eff (TDCF - 1) 

7dcf^^ sin/ 



(19) 
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values) than the initial pericentre. We found that initially retrograde orbits 
became captured in this model. This feature can be naturally explained as 
follows. Dissipation of the orbital energy in each intersection with the disc 
tends to increase the binding energy of the orbiting object. In the case of 
originally prograde orbits, however, the orbiter obtains sufficient amount of 
angular momentum, which saves it from being captured by the black hole. 
Finally, the object settles into a circular orbit in the disc plane. However, an 
object which started with retrograde orbit does not acquire enough angular 
momentum during the period of nearly perpendicular intersections with the 
disc. Due to continuous losses of energy it typically becomes captured by the 
hole. 

To clarify previous results based on the relativistic treatment we compared 
them with the corresponding Newtonian "elliptic" mapping (see Appendix). 
To be consistent, we also reduced formula for the star-disc interaction by 
eliminating the Lorentz factors 7dcf in (19) and instead of the Lorentz boost 
from LNRF to DCF we used the Galilean transformation. Fig. 3a shows 
the fully relativistic model with the Schwarzschild background metric, while 
Fig. 3b is the Newtonian analogue. We have chosen formally the same initial 
eccentricity and inclination in both Figures: e ~ 0.83 and / = 130°. The orbit 
is initially retrograde and in the relativistic case it becomes captured by the 
central black hole. On the contrary, this docs not occur in the Newtonian 
case and the orbit is circularized to some definite radius. (In a realistic case, 
however, the orbiter can be tidally disrupted before it is captured but this 
depends on its internal structure and we do not consider such possibility in 
the present paper.) Because the interaction of the orbiter with the disc was 
always chosen to be weak, the time scale for the precession of the pericentre 
is much shorter than the time scale for the evolution of the other orbital 
parameters. As a consequence, the points of intersection with the disc fill 
the interval between the current pericentre and apocentre in the relativistic 
case (Fig. 3a). There is no precession of the pericentre in the Newtonian 
case and thus only an 'adiabatic' evolution of the orbital parameters is seen 
(Fig. 3b). 

We also show the results of integration with the nearly extreme Kerr 
black hole. In this example we have chosen a = 0.9981 for dcfiniteness 
(Thorne 1974), the same initial eccentricity and inclination, and the same 
initial apocentre and pericentre (expressed in the gravitational radii) as in 
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Fig 2. We observed (Fig. 4) a slightly shorter circularization time than in 
the corresponding Schwarzschild case, but the qualitative features of the disc- 
orbit interaction remained unchanged. They may be changed, however, when 
details of the structure of the accretion disc are taken into account because 
the structure of the disc and the location of its inner edge depend significantly 
on a. Again, we found that initially retrograde orbits get captured by the 
black hole. 

4 CONCLUSIONS 

We assumed that the low-mass compact object interacts with the thin accre- 
tion disc twice per each revolution — exactly when it crosses the equatorial 
plane of the black hole (impulsive approximation). We described the relevant 
equations and we employed them in the fast numerical code computing the 
evolution of the trajectory. We found that the effective time of circularization 
is shorter than the time to change the orbital plane into the plane of the disc. 
This conclusion is in accordance with previous results of Syer et al. (1991) 
based on Newtonian gravity. However, we have also observed short periods 
during the evolution when eccentricity increases. In particular, this increase 
occurred in model described by Eq. (18) during the transition period where 
the initially retrograde character of the orbit is changed to a prograde one. 

The star-disc interaction was described by a phenomenological parame- 
ter characterizing the magnitude of the change of orbital parameters in each 
collision. This phenomenological description is satisfactory provided the disc 
remains thin and the orbital parameters are changed only at the moment of 
transition of the star through the equatorial plane of the central black hole. 
A description of the interaction which would result from a detailed physical 
model is not crucial in this case; we want to improve our understanding of the 
interaction in future work. Effects of the dynamical friction and direct accre- 
tion acting on an object moving through the gaseous medium were studied 
by a number of authors under various conditions (recently by Petrich et al. 
1989). In our highly supersonic and turbulent case, the approach outlined 
by Zurek et al. (1992) appears the most appropriate one. 

One can specify parameters of the model for the case of NGC 6814. Each 
single long-duration observation of EXOSAT or Ginga covers less than 30 
revolutions of the orbiter. The characteristic time-scale for the precession 
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of nodes is much longer than the orbital one. The estimate which adopts 
the maximum value of a = 1 and radius of the orbit ~ 50 gravitational 
radii of the central black hole yields the ratio of the Lense-Thirring to the 
orbital frequency ~ 0.005 . This means that the orientation of the orbit is not 
significantly changed during each observation. However, the interval between 
EXOSAT ohsBTvaiions and Ginga observations was certainly long enough and 
a resulting change in the orientation suggests a possibility to understand the 
perfect stabihty of the orbital period detected by both satellites and at the 
same time the puzzling change in the light curve profile. We do not want 
to speculate further on this important subject until the star-disc collisions 
and dynamical friction acting on the star are better understood (work in 
preparation). 

To conclude, adopting the model of star-disc interactions as an explanation 
of the origin of the NGC 6814 hght curve, we see one important contribution 
from general relativistic effects which is due to the pericentre shift and Lense- 
Thirring precession. These effects drag the point on the orbit where the star 
crashes through the disc. They also modify the velocity at which the star hits 
the disc as well as the orientation of the orbit with respect to the observer. 
This fact has two consequences: (i) additional periodicities corresponding 
to the precession frequencies are present and can potentially be revealed in 
the power spectrum of the signal from the source (Karas & Vokrouhlicky 
1993); (ii) long-term evolution of marginally stable and marginally bound 
orbits is very different compared to the orbits with identical initial parameters 
treated in Newtonian theory of gravity. The first consequence above gives 
us a possibility to detect Lense-Thirring precession induced near the core 
of an AGN. If the corresponding frequency is not present we will be able 
to conclude that the central black hole (if any) is non-rotating. This would 
be extremely important information especially from the point of view of 
electromagnetic scenarios of AGN which require a rotating black hole. The 
second consequence is particularly important in describing the capture of the 
star into a bound orbit around the central black hole. Although many models 
assume a star located on such an orbit, the very process of the capture is not 
well understood. 
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APPENDIX A: THE MAPPING ALGORITHM 



This Appendix outlines the exphcit form of the mapping (r, 0; ri)i — * (r, 0; r]) f 
from Section 2. Ahhough our approach is straightforward, we beheve that it 
has not yet been employed by other authors. As we found it very advanta- 
geous for practical purposes, we describe derivations relevant for this work in 
some detail. We constructed a numerical code which employs efficient rou- 
tines for evaluation of elliptic integrals and Jacobian elliptic functions (Press 
et al. 1986, Press & Teukolsky 1990). The code achieves better precision 
and about two orders of magnitude higher speed compared to direct numer- 
ical integration of the geodesic equation in its equivalent form of first order 
differential equations. (We used direct integration to check the code.) The 
two cases, < a < 1 and a = 0, are technically somewhat different and we 
discuss them separately. 

Al The case < a < 1 

(i) Evaluate the latitudinal integral between two successive intersections with 
the equatorial plane: 

where K(k) denotes the complete elliptic integral of the first kind. 

(ii) Distinguish the three cases which may occur: 

• Case I — four real roots of R{r) = 0, r4 < r < r^. 

• Case II — four real roots of R{r) = 0, r2 < r < ri. 

• Case III — two real and two complex roots R{r) = 0, r2 < r < ri. 

Denote r^j = 1 (r^j = —1) if r is increasing (decreasing) at rj; analogously r]f 
for r/. 

(iii) Evaluate the increase of between each two radial turning points: 

KK{ki) (Case I and II), 

6Ir = I (A2) 
, ^ K(k2) (Casein), 
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where 



{ri - r2)(r3 - r4) 
{ri - r3)(r2 - ' 



k2 

= (xi - n)^ + xl , 

K 



Apq 



(iv) Denote 



K,F{(p, ki) (Case I and II), 

F{^p,k2) (Casein), 

if 77/77j > 0, 



' m5Ir- Tjff 



Ir — i 



_ (m - 77/) 5/^ + r]ff if 77/77j < 0. 



Here, 



sin^ if 



( (n-r3)(r.-r4) (Casel), 



(ri-r3)(r-,-r-2) (Qasell), 



(Case III); 



(A3) 



2 g(ri - Ti) 

F{ip^ k) is the incomplete elhptic integral of the first kind and m is the number 
of turning points in r between the two successive intersections. In addition, 
one has to check whether the trajectory still remains above the horizon if the 
lower turning point is located below the horizon. 

(v) The radial coordinate of the intersection is 



with 



r irr^r,)r^(r,-r,rna ^QaSel), 
ri-r3+{r3-r4ja V /' 

(ri-r3)r2 + (ri-r2)r3g fQagg J J) 



a = sia^{u, ki) , 



(A4) 
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u 



or 



with 



qria + pr2 
p + qa 



(Case III), 



(A5) 



[1 + cn{u, k2)] 



sn{u, k) and cn(M, k) are Jacobian elliptic functions. At this point we are 
able to compute the r-coordinates of the intersections, which is sufficient to 
determine the evolution of eccentricity and inclination of the orbit and the 
number of revolutions before the trajectory becomes captured by the black 
hole or escapes to £^ > 1 (and then presumably to infinity); (j) coordinates are 
also needed if we wish to study the precession. Finally, we need coordinate 
time to relate the revolutions to time as measured by a distant observer. In 
the Case III the orbit is in practice captured by the black hole after a few 
revolutions. Thus we exclude this case from further considerations. 

(vi) Evaluate the following quantities. 

Case I: 



I± = [{r4 - ri)n {ip, n±, fci) + (r± - r4)F{ip, ki)] + /±, (A6) 

(ra -r4)(r± - ri) 



(ri - ra) (r± - r4) 



Jr + 2Kr — Kr4 



— o?- 



2r4ai- 



,a —a 



a" 



+ Jr + 2Kr, 



(A7) 



with 



U = F{^,k^), 
Vi^U (ifi, -a\ ki) , 
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Vo 



2(a2_i)(A;2_«2; 

kl -a^)U- sn([/, k,) cn([/, k,) dn([/, k,) „ . , 

^ 1 — a^SB.^{U, kl) 



a 



K± 



ri - ra ' r4(ri - rs) ' 

sm = , 

(ra - ri)[ri - Ti) 

2 

(r± - ri)(r± - ri)^{ri - rz){r2 - r4)(l - 



Case II: 



/± = [(r2 - r3)n ((^, n±, /ci) + (r± - r2)F((^, /ci)] + /±, (A8) 

(r2 -ri)(r± - rg) 



(ri - ra) (r± - r2) ' 



Jr + 2^^ = >^f2 



l4 + 2'A\u + 2^^V,+ 



,a — a 



+ Jr + 2Kr 



with U, Vi and V2 defined as above, and 



n - r2 
n - ra 



rsjn - r2) 
^2(^1 - ra) ' 



• 2 ('^i - ^3)(rj - r2) 

sm ^ = r , 

(ri - r2)(ri - ra) 

2 



(r± -r3)(r2 -r±)W(ri - r3)(r2 - r4)(l - i:^) 



(A9) 



E{ip^ k) and 11 ((/?, n, k) are incomplete elliptic integrals of the second and the 
third kind, respectively. Integration constants /±, Jr and are the values 
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of integrals in (15) and (17) evaluated between and the last turning point. 
Thus, they depend on the number of turning points in r (m) and the sign 
of initial radial velocity (rji) and they can be given in terms analogous to 
Eq. (A3) . We skip explicit expressions because they are rather lengthy. 

(vii) Finally, 

J, = , , (AlO) 

= H c2 - (All) 

ay 1 — ^ 

with E{k) and n in. k) being complete elliptic integrals of the second and 
the third kind. Now we have all necessary quantities required for complete 
mapping of relevant trajectories in the Kerr metric. 

A2 The Schwarzschild case (a = 0) 

We present the case of the Schwarzschild background metric separately, even 
though the general formalism developed for the Kerr metric can also be 
applied. The reason is twofold: (i) the formulae valid for the general Kerr 
metric often include the angular momentum parameter a in denominators 
[e.g. Eqs. (AlO), (All)]. These apparent singularities cancel out in the 
limit a — > 0, but they are the source of difficulties in numerical evaluation; 
(ii) as the symmetry of the space-time is now higher and geodesies in the 
Schwarzschild metric remain always planar, we can avoid integration of the 
latitude 9, restricting ourselves to the current orbital plane of the test particle 
spanning one loop of the trajectory above/below the disc (the true orbital 
plane of the object is changed due to the interactions with the disc). Thus 
we reduce the order of the mapping by evaluating the integrals in polar 
coordinates in the orbital plane. The orbital plane differs from the disc plane 
by the current value of the inclination. 

We consider the following coordinate systems: (i) Schwarzschild spherical 
coordinates {r,9,(f)); latitude 9 is measured from the axis of the disc plane 
and polar angle (j) in the disc plane {(j) — direction can be chosen arbitrar- 
ily), (ii) (r,'&) polar coordinates in the current orbital plane of the orbiting 
object where the angle is measured from the actual nearest preceding apoc- 
entre of the unperturbed trajectory with the current orbital parameters. Let 
us clarify better the concept of the origin, as it is intimately connected 
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with our technique. The analytic integration of the geodesic motion in the 
Schwarzschild space-time is advantageously carried out if the polar angle 
in the orbital plane is measured from the nearest preceding apocentre. In 
each step of the mapping procedure we are interested only in one orbital loop 
above/below the accretion disc; then the interaction with the disc changes the 
orbital parameters for the successive loop. It is this orbital loop, where the 
free motion of the test particle in the Schwarzschild background is applied. 
However, the orbital loop which is under consideration may not necessarily 
contain the apocentre of the orbit. Thus, apart from the true trajectory of the 
object we introduce a reference trajectory of the object with the same orbital 
parameters as for the true one and coinciding with the true trajectory just on 
the current segment. This fictitious reference orbit defines the d origin — it is 
measured from the nearest preceding apocentre of the reference trajectory. 

The equations of motion covering a single mapping step are as follows 
(e.g. Chandrasekhar 1983): 



and u = 1/r. Constants of motion are defined as £ = —pt and C = [note, 
that angular momentum C is defined with respect to the fictitious orbital 
plane, not with respect to the disc plane hke $ in the Kerr case]. They are 
related to the tetrad components of the four-momentum in the locally static 
frames by means of Eqs. (6)-(7). We define / = | — /9 as the inclination of 
the fictitious orbit with respect to the fixed reference disc plane [it can be 
equivalently expressed using the LNRF tetrad components in the equatorial 
plane: tan7 = p7P^; cf. Eq. (8)]. 

Again, we will concentrate on orbits characterized hy £ < 1 and £ 7^ 0. 
The signs of the first and the last term of the polynomial U (u) guarantee 
at least one positive root of the equation U{u) = and, as U{u = 0) < 0, 
we conclude that this root corresponds to the apocentre of the orbit. The 
type of the orbit is determined by the properties of the other two roots of 




(A12) 




(A13) 



where 



U{u) = 2u^ -u' + 2C-^u - (1 - £^)C-\ 
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the equation U{u) = 0. The roots cannot be real and negative at the same 
time (Chandrasekhar 1983). We exclude the possibility of multiple roots of 
as before. Hence, we are left with the two kinds of orbits characterized by: 

• the three positive real roots of U{u) = which we arrange according 
to magnitude: ui < U2 < Ms, 

• one positive real root (tii) and two complex conjugated roots {uc,v^) 
of U{u) = 0. 

The first item still encompasses two types of orbits: (i) those captured by 
the black hole in the sense that they have no pericentre above the horizon 
(the apocentre of such orbits is always less than 6 — where the last stable 
orbit and presumably the inner edge of the accretion disc are located in 
the Schwarzschild case; hence we do not consider these orbits); (ii) quasi- 
elliptic orbits bound between the turning points, Ui and U2; we will call them 
Case I orbits. In terminology used by Chandrasekhar (1983) our Case I orbits 
correspond to the orbits of the first kind. We call the Case II orbits those of 
the second item above. They have no pericentre and fall unavoidably to the 
black hole (they correspond to the orbits with purely imaginary eccentricity 
in Chandrasekhar 's terminology). 

One can easily find a simple rule dividing both cases of the orbits: < 12 
implies the Case II orbit. The Case I orbits are characterized by > 12 
and simultaneously 

{l-2u.){l + C\l)>£', 

where = l + \C\~^\/£,'^ — 12. In what follows, we will describe an algorithm 
for the mapping of these two cases of orbits in detail. We will pay special 
care to the Case I orbits, as they will be shown to be the most important in 
astrophysical applications. 

A2.1 Case I orbits 

In close analogy to the Newtonian case, Eq. (A12) is advantageously inte- 
grated in terms of the relativistic "true anomaly" 

u{x) = /^(l + ecosx), 

where 

Ui +U2 U2- Ui 

ji — — , e = . 

2 U2 + U1 
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The quantity e can be interpreted as the eccentricity of the orbit. We do 
not write the exphcit form of primitive functions obtained by integration 
(Chandrasekhar 1983) but we give formulae for the mapping which we need 
in our present work. After some manipulation we arrive at the following form 
of the mapping: 



Uf ^ U2- {U2 - Ui] 

where now 



1 - A;2$$^ 



(A14) 



Ui 



Ui 



U2 - Ui 



sn 



TT 



a;, k 



U2 



Ui 



U2 - Ui 



cn 



TT 



cu, k 



UJ 



— \/Y^^2u2--'^u\^ — 2{u2 — ui)u~ 



Mapping of the sign function rj is given as follows: 



Vf 



sgn[a{xi) + - a{0)] if r^i = -1, 



where 



Xi 



sgn[o-(xi) - tt] 

= arccos($ — \E') 
a{xi)^2u;-'F 



TT 



if Vi = 1, 

(0,7r), 

- Xi 



(A15) 



'-,k 



It is instructive to discuss the Newtonian limit of the mapping formula 
(A14) in which the terms proportional to some power of 1/c are neglected. 
This limit is now obscured by the fact that wc imposed widely used "relativis- 
tic convention c = 1" while now we want to suppress the terms containing 
c in the denominator. Careful bookkeeping of c in the preceding equations 
suggests that the Newtonian hmit corresponds to the fictitious procedure: 
a;— >0. As a result we obtain 



U 



N 
2 ' 



Uf = —Ui + U-^ 

where the Newtonian boundaries are found to be 

(l ± Vl + 2E^C^) , 



(A16) 



u 



.N 
1,2 
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(we retain the negative value for the Newtonian energy of bound orbits, as 
seen from the above definition of S'^). Eq. (A16) is the correct expression 
for the "eUiptic" mapping. Formula (A16) is surprisingly simple, showing 
the linearity of the elliptic mapping. Moreover, we also have the simple rule 
rjf = —rji. By contrast, the full relativistic mapping (A14) for the orbits of 
the Case I is highly nonlinear and the mapping in u coordinate is coupled 
with the mapping of 77 function due to r]i in Eq. (A14) and u in Eq. (A15). 

Let us turn to Eq. (A 13) describing the mapping in the t-coordinate. We 
start by expressing the indefinite integral on the right-hand-side. Changing 
the ■& variable to x according to the relation 

^ = -uj^l-k^coB\xm ^ -a;A(x, k) , 
one arrives at the primitive function 

I ^ n ( A; 

[l-2/.(l-e)] V 2 '2Ml-e)-l' \ 



1 

4 



2/^2(1 -e2)(l + eA;2) 



A(x, fc) sinx 
2[e + cos2(x/2)] 



+[1 + 2e(l + k"") + 3e'A;2]e-^n (^^^, k 



(A17) 



where 



Ui 

e — 



U2 - Ui 

Now, the algorithm for evaluating the time step associated with the mapping 
is as follows: 

rii = —1: introducing Q — sgn{a{0) — cr{xi) — we obtain 

f.-f - / ^(^/) - ^(^^) if C.>0, . . 

2T(0) - T(x/) - T{x^) if < 0, ^""'^^ 
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where X/ is determined by relations 




1^ 



Tixi)-T{xf) if Ci>0, 
T{xi)+T{xf) if Ci<0, 



— 1: introducing Q = sgn(a{xi) ~ ^r) we obtain 
tf -ti = 

and now X/ is determined by relations 



(A19) 



f sin 




Cisn 










[ 


[ cos 


I J 





The Newtonian limit of Eq. (A 18) is 



tf - ti 



1 



arctan 



el sin | 



+ 



e sin 1?,: 



1 — cos^ 



(A20) 



where 



cosi?i = $-*, sinT?j ^ rjiy/l - - . 
and functions $ and are defined as in (A14) and evaluated for ui 

N 

U2 = U2 ■ 



= U 



N 
1 ' 



A2.2 Case II orbits 

These orbits have no pericentre. They are captured by the black hole in most 
cases, even though the interaction with the accretion disc can in principle 
modify the orbital parameters and change the type of the orbit. Conse- 
quently, we do not perform the mapping in full detail — wc skip the deriva- 
tion of the time interval {tf — ti) which will not be needed for this type of 
orbits. We still need to describe the complete algorithm for the mapping in 
the w-coordinate and associated r^-function. We were unsuccessful in finding 
a compact expression for the mapping Ui — > Uf{ui]Ui,Uc) similar to that 
presented in formula (A14), and thus in the following we give the algorithm 
of the mapping in several steps (well suited for programming) . 
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(i) Evaluate the following quantities: 



/x=-(1-2mi), e = ^3 - /(I - , 



+ 2mi - 1 1 



25 ' 



e(l - 2mi1 



1 — if Ui> ji, 
\J^- if Ui < II, 
and the angle Ui e (— f , f ) which is defined by 

sin^ = 1 - 27~^ec pe/xe^ + (6/^ - IjeJ , 

where sgn(a;j) = sgn(es — . 

(ii) Exclude the captured trajectories by evaluating 

t?, = [ir(fc+)-F(a;,,fc+)]rV2; 

the particle will be captured by the black hole before reaching the equatorial 
plane if < tt and rji — —1. Otherwise we continue to the following step. 

(iii) Define 



and 



A* = F{u, k+) - r],ii5^'\ = sgn(i^(A;+) + A.) , 



CiSn(A*,/c+) 
cn(A^, k+) 



\ Us 


1=1 


I '^^ J 





(iv) Finally, the mapping which we look for is given by 

Uf^n(l + e7[sg„(«,_5)i) , Vf = CiVi, 
where we have used 

2e/x ± •^4e2/i2 + ^^^^ {^+kI - 7_) 



(A21) 



1..-2 

c 
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5 = -^l-2(6/i-l)7;^ 

These orbits are unstable in the sense that they reach the singularity at r = 
and thus in the Newtonian limit there are no analogous orbits corresponding 
to this case. 

Eqs. (A14) and (A21) give the mapping for astrophysically interesting 
cases of geodesies in the Schwarzschild geometry. 

FIGURE CAPTIONS 

Figure 1. The sequence of successive intersections with the disc. We plot 
radial coordinate r of the intersection on the ordinate and the number of 
intersections N on the abscissa. For > 5x 10*^ the upper and lower bound- 
aries of the distribution of intersections in the Figure get closer to each other, 
which means that the eccentricity of the orbit decreases; simultaneously as 
the orbit is circularized it is also ground to the plane of the disc. In this 
case, (X = and = 0.9999. The initial pericentre distance is 7, eccentricity 
0.7, inclination / = 103°. The arrow indicates the moment when the orbit 
changes its character from retrograde to prograde (/ = 90°). 

Figure 2. As in Fig. 1 but for the second model of the star-disc interaction 
[Eq. (40)]. Initial pericentre distance is 30 and eccentricity e = 0.83. Two 
initial inclinations are compared — (a) / = 35°, (b) / = 80°. The proportion- 
ality constant in (40) is taken 10~^. Originally less inchned orbit (a) settles 
on the circular orbit in the disc with radius of ~ 53, while the more inclined 
orbit (b) grinds to the circular orbit at radius ~ 17.7. We verified, as an 
example, that keeping the initial pericentre at 30 the results are not very 
sensitive to the initial eccentricity provided it is ^ 0.75. 

Figure 3. The graph (a) shows a similar orbit as arc those in Fig. 2, but now 
the initial trajectory is retrograde with an inclination of 130°. As commented 
in the text, it is captured by the central black hole. The graph (b) shows the 
'Newtonian analogue' of (a). Only the 'adiabatic evolution' of the apocentre 
(upper curve) and the pericentre (lower curve) is seen and there is no shift 
of pericentre. 

Figure 4. Orbits analogous to those in Fig. 2, but with a nearly extreme 
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Kerr metric (a = 0.9981). The initial pericentre distance is 15 (the same as 
for the Fig. 2 orbits if expressed in gravitational radii of the central black 
hole). The initial inclinations are again chosen to be / = 35° (a) and / = 80° 
(b). 
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Relativistic precession of the orbit of a star 
near a supermassive black hole 



ABSTRACT 

Wc study the gravitomagnetic (Lense-Thirring) precession of the trajectory 
(approximated by a geodesic) of a star orbiting a supermassive rotating 
(Kerr) black hole. We do not assume any particular value for the eccentricity 
or inclination of the orbit or the angular momentum of the black hole. We 
also discuss the periodicity related to the relativistic shift of the pericenter. 

The Seyfert galaxy NGC 6814 is an example of the object for which effects of 
relativistic precession could be detectable and we discuss the relevant preces- 
sion frequency for this case. The remarkably stable phase of several patterns 
and their position in the light curve impose strong restrictions on the model 
of this object. We conclude that, according to our present knowledge, a star 
colliding with an accretion disk is somewhat improbable though not com- 
pletely excluded as a model of NGC 6814. Our arguments are independent 
of studies of the period stability. 

Subject headings: galaxies: active — galaxies: individual (NGC 6814) — 
general relativity — black holes 

1. INTRODUCTION AND MOTIVATION 

It is widely accepted that supermassive black holes (SBH) are located in cores 
of active galaxies and quasars (Begelman, Blandford & Rees 1984; Shlosman, 
Begelman & Frank 1990). The mass of the SBH is usually estimated to be 
in the range M ^ 10^-10^^ Mq. Anomalous energy output of active galactic 
nuclei (AGN) may result from an accretion process: the matter is attracted 
from the surroundings of AGN or it comes from tidally disrupted stars passing 
too close to the SBH. The accretion disk is formed and the matter eventually 
falls onto the black hole. Although this scenario can be called standard, the 
evidence for both accretion disks and SBH in AGN is only indirect (Frank, 
King & Rainc 1985; Blandford, Nctzcr & Woltjcr 1990; Falckc et al 1993). 
The best evidence comes from studies of the central surface brightness of 
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the nuclei, stellar velocity dispersion, spatial distribution of stars and X-ray 
emission (Young et al. 1978; Sargent et al. 1978; Binney & Tremaine 1987; 
Dressier & Richstone 1988, 1990; Kormendy 1988a, b; Halpern & Filippenko 
1988; Dressier 1989). General relativistic effects may have important conse- 
quences for the axial symmetry and stability of accretion disks (Bardeen & 
Petterson 1975; Abramowicz 1987). However, the presence of black holes 
in AGN is largely masked by violent plasma processes in the surrounding 
medium. Electromagnetic models of energy extraction assume that the SBH 
rotates (Blandford k Znajek 1977; Macdonald & Thorne 1982; Phinney 1983; 
Kaburaki & Okamoto 1991; Okamoto & Kaburaki 1991 and references cited 
therein). In this scenario, the energy of an AGN comes, at least partially, 
from the rotational energy of the SBH. Evolution of the angular momen- 
tum of the black hole under such process was studied by Park & Vishniac 

(1991) . It appears that the angular momentum determines the energy output 
of the AGN in a nontrivial manner, with the maximum at some particular 
value (Bicak & Janis 1985). Unfortunately, it is unclear how the value of 
the angular momentum of the SBH could be determined by an independent 
observation. The present paper deals with this problem. 

It has been proposed that a star could be captured in a bound orbit around 
the SBH by the tidal distortion and associated dissipation of energy (Fabian, 
Pringle & Rees 1975; Frank & Rees 1976; Press & Teukolsky 1977; Lee & 
Ostriker 1986; Rees 1988), tidal disruption of a binary star (Hills 1988) or a 
cluster (Novikov, Pethick & Polnarev 1992), or by cumulative effects of inter- 
actions with an accretion disk (Syer, Clarke & Rees 1991). Various aspects 
of star-disk collisions were studied by Ostriker (1983), Zentsova (1983), Syer 
et al. (1991), Zurek, Siemiginowska & Colgate (1992), Sikora & Begelman 

(1992) and Vokrouhhcky & Karas (1993). We assumed that the SBH forms 
such a binary system with a low mass star in an orbit inclined with respect 
to the plane of an accretion disk. The disk is presumably thin and its axis is 
aligned with the rotation axis of the black hole (Bardeen & Petterson 1975; 
Kumar & Pringle 1985) but the model can easily be generalized to a more 
complicated geometry. Radiation from the disk is periodically modulated 
each time the star crosses the disk. We did not attempt to specify any par- 
ticular mechanism of the modulation. We also ignored secular changes of 
the orbital parameters of the star due to collisions with the disk and cor- 
responding contributions of the collisions to the precession frequency. We 
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concentrated on a gravitomagnetically induced precession of the orbit of the 
star. This precession is a general relativistic effect; it occurs if the central 
black hole rotates and it becomes important when the star revolves close to 
it. In particular we attempted to pick up relevant frequencies in the power 
spectrum of the simulated signal which are independent of very complex and 
poorly understood details of the star-disk interaction. Given the model of 
interaction our approach can easily be adapted. We took into account all the 
relativistic effects affecting the motion and energy of photons arriving from 
the source to a distant observer. 

We assume that the star moves along a geodesic around a Kerr black hole. 
The gravitomagnetic effect was originally treated by Lense & Thirring (1918) 
in the weak-field limit and by Wilkins (1972) in the case of a spherical or- 
bit, r = const, around a black hole with the extreme value of the angular 
momentum parameter, a = 1. According to our knowledge, a generalization 
to an eccentric and inclined orbit around a black hole with arbitrary value 
of parameter a e (0, 1) has not yet been discussed in the literature. As a 
consequence of the gravitomagnetic effect, orbital nodes are dragged in the 
sense of rotation of the black hole. This dragging affects the position of the 
source with respect to a distant observer. Wc try to extract the information 
which might help us to determine whether the central SBH rotates or not — 
provided other details of star-disk collisions are known. The weak-field limit 
of the angular velocity of the gravitomagnetic dragging is 

fiLT-4xlO=§^ s-', (1) 

where r and a are the radius of the orbit and the angular momentum pa- 
rameter of SBH measured in dimensionless geometrized units, r = 6.7 x 
10~^(f/l cm). Several experiments to confirm the dragging effect in the limit 
of a weak gravitational field have been proposed but they have not been 
carried out as yet (Everitt 1974; Will 1981; Braginskij, Polnarcv & Thorne 
1984; Ciufolini 1986). The strong field regime of the spin-orbital interaction 
is also intensively investigated for a relatively broad class of gravity theories 
within the framework of the parametrised post-Keplerian formalism. The 
main applications of this approach are directed at the interpretation of the 
binary pulsar data (e.g. Damour & Taylor 1992). Though not yet con- 
firmed by direct observations, the gravitomagnetic effect is considered as a 



32 



firm consequence of general relativity. In the strong-field region close to the 
black hole the gravitomagnetic precession depends on four parameters — the 
pericenter distance of the orbit, eccentricity of the orbit, inclination with 
respect to the equatorial plane of the black hole, and, of course, the angular 
momentum parameter. The above mentioned processes of tidal interaction 
can set the star on an initially very eccentric trajectory and we therefore 
could not restrict ourselves to circular orbits. Star-disk collisions (Syer et al. 
1991; Vokrouhlicky & Karas 1993), tidal effects (Press, Wiita & Smarr 1975; 
Lecar, Wheeler & McKee 1976; Boyle & Walker 1986; Zahn 1977 and 1989; 
Zahn & Bouchet 1989; Tassoul 1988) and gravitational radiation (Peters & 
Mathews 1963; Zeldovich & Novikov 1971) tend to gradually circularize the 
elliptic orbit, however. Naturally, these effects modify the precession fre- 
quency, but we assumed that the star is a dwarf or a low-mass compact 
object and ignored the corrections. We also ignore all possible effects of 
the magnetic field on the space-time geometry. For the discussion of exact 
solutions of Einstein-Maxwell equations describing a black hole in a mag- 
netic field cf. Ernst (1976), Karas (1991), Manko & Sibgatullin (1992), and 
references cited therein. 

In other words, we ignored all effects which could induce the precession of the 
star's orbit except the gravitomagnetic effect. We believe it is an appropriate 
approach in solving the problem if these effects (like star-disk interactions or 
gravitational radiation) are also weak. We will discuss other effects elsewhere 
so that the whole subject is treated in steps. 

We determined the value of the precession frequency (or alternatively the 
nodal shift per one revolution, 5(f)) by direct integration of the geodesic 
equation in terms of elliptic integrals. Gravitational radiation losses can 
be neglected provided the change of the energy is small; for a circular orbit 
with energy £ one obtains the dimensionless estimate (e.g. Rees, Ruffini & 
Wheeler 1974) 

I 6.4 - (-}] «1, (2) 



where is the mass of the star. Suppose we estimate M, r and an upper 
limit on the change of orbital period P from observation. Then equation 
(0) provides us with the upper limit on the mass of the companion star [cf. 
Sikora & Begelman (1992) and King & Done (1993) for the application to 
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NGC 6814]. Analogously, if R^, is the radius of the star one can obtain a 
constraint on the pericenter distance Rp which is based on the tidal limit 
(Carter & Luminet 1983; Luminet & Marck 1985; Rees 1988): 



This paper is organized as follows. In Section 2 we derive the azimuthal 
shift of orbital nodes. An unambiguous value can only be given in the case 
of spherical trajectories. The shift oscillates between the maximum and the 
minimum values, ^0max and 50min5 if the orbit is eccentric. We define a suit- 
able probability distribution for 50 and compute, numerically, a mean value 
(50) which determines the gravitomagnetic precession averaged over a large 
number of revolutions of the star on an eccentric trajectory around the SBH. 
Corresponding formulas are rather complex and, therefore, we also present 
a table of numerical values in Appendix. In Section 3 we present simple 
examples to illustrate how the gravitomagnetic frequency could be extracted 
from observational data. (In the present contribution, we use simulated data 
from a simple model rather than real data.) Naturally, the frequency cannot 
be too small so that the data cover at least several periods if the effect is to 
be detectable. The estimate of the time interval between successive collisions 
with the disk [see eq. (pO|) below] leads us to assume that the star crosses 
the accretion disk near the horizon (typically a few tens of the gravitational 
radius of the SBH) in the region where the collisions can modulate the disk 
radiation in the optical/X-ray bands. Possible physical mechanisms for the 
modulation were discussed by Mushotzky (1982), Guilbert, Fabian & Ross 
(1982) and Zentsova (1985). The flux of radiation is modulated by the orbital 
period (short time-scale) and by the perihelion shift and the Lense-Thirring 
period (longer time-scale). In order to compute observable effects we consider 
both the time delay and the focusing of photons coming from the disk region 
to a distant observer. The focusing effect strongly enhances the influence of 
the Lense-Thirring precession for observers with large inclination (edge on, 
with respect to the accretion disk). Finally, we apply the model to the case 
of the Seyfert galaxy NGC 6814. Periodic modulation of the X-ray emission 
from this galaxy has been confirmed on the time-scale ~ 12, 200 s (Mittaz & 
Branduardi-Raymont 1989; Fiore, Massaro & Barone 1992; Done et al. 1992) 
and possible mechanisms were proposed by several authors (Abramowicz et 
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al. 1989 and 1992; Done et al. 1991; Honma et al. 1991; Wallindcr 1991 and 
1992; Abramowicz 1992; Bao 1992a, b; Sikora & Begelman 1992; Rees 1993; 
Vio et al. 1993; King & Done 1993). 



2. PRECESSION FREQUENCY— DETAILS OF THE 

CALCULATION 

2.1 Gravitomagnetic precession 

Geodesic motion in the Kerr space-time can be integrated in terms of elliptic 

integrals (Carter 1968). The appropriate form can be found in Vokrouh- 
hcky & Karas (1993). We use the standard notation for the Kerr metric 
(Bardeen 1973) and we refer the reader not familiar with the Kerr metric to 
this reference. We assume that the motion of the star is quasi-elliptic with 
the pericenter outside the black hole horizon, r = = 1 -|- y/1 — a^, and 
with positive energy, < S < 1. The locations of the pericenter r = Rp and 
the apocenter r — coincide with the two upper roots of the polynomial 

R{r) = - l)r^ + 2r^ + [(£' - l)a' - - Q] + 2/Cr - a'Q, (4) 

where £ = —pt, $ = and IC are the usual constants of motion, Q = 
JC + — aSy. In our case, all roots of R{r) are real. We denote them by 
R3. > R^ > Rd, ^ Ra'i R?, > R+- Other possible combinations — e.g. two 
complex conjugated roots — are exchided by our previous assumptions on 
energy and location of the pericenter above the horizon (Stewart & Walker 
1973). We further assume that the star periodically crosses the equatorial 
plane, 6 — 7r/2; latitudal motion is restricted by | cos 6' |< //■-, where is 
the lower of the two positive roots ix± of the polynomial 

Q{^x) = (1 - / - [Q + a" (l - + /.^ + Q; (5) 

For r ^ one can interpret arccos as the inclination of stellar orbit. 

A non-equatorial geodesic around a rotating black hole is not planar; inter- 
sections with the 9 — 7r/2 plane are dragged in the sense of rotation. At first, 
we apply the mapping 

[r, 0, t, sign(r)]^ — > [r, 0, t, sign(r)]^^^ (6) 
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which expresses coordinates of the {n + l)-st intersection in terms of coor- 
dinates of the previous, n-th intersection with the equatorial plane. The 
azimuthal shift of the orbital node 5(f) per one revolution is then defined by 

S4> = (f)n+2 -4>„-2tT (7) 

(0-coordinate is not restricted to (0, 27r) in this convention). One can find 

^ _ (-Ra — R3)Rp + (Ra. — Rp)Rzcr 

Ra — R3 — [Ra. — Rp)<y 

0n+i-0n = [2{aS - + /++[2(a£ - + /_+$J^ , (9) 

tn+1 -tn = £{Jr + 2Kr) + 2 [B+R+S + a(a^ - ^)A+] 1+ 

+2 [B_R_S + a{aS - ^)A_] I_ + ASI^ + a^£K^ . (10) 



We denoted 



(7 = sn^(ii, /ci) , (11) 



1 

^=2 



(i?a - i?3)(i?p - Ra){1 - ^2) {l, - ir) , (12) 

(i?a-i?p)(/?3-i?4) 



{Rn — R3)iRp — R-i) 



R-^- — R— R-\- — R— 

R± ^i±VT^, 

l=^^(W/^+), (13) 

I± = «± [{Rp - i?3)n ((/^, n±, A;i) + (i?± - Rp)F{ip, h)] + /±, (14) 

^ _ (-Rp — Ra){R± — Rs) . 
(i?a — -Ra) (R± — Rp) ' 

F((/7, A;) is the incomplete eUiptic integral of the first kind, K{k) = F(7r/2, k) 
is the corresponding complete integral and sn{u, k) is the Jacobian elliptic 
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function (e.g. Byrd & Friedman 1971; Gradshteyn & Ryzhik 1980). Anal- 
ogously, n((y9, n, k) is the elliptic integral of the third kind. In equations 



a'* 



+ 2^ f/ + 2 



c? — c? 



Vo 



+ Jr + 2Kr 



(15) 



with U, V\ and V2 defined by 

f/ = F(<^,fci), 

Vi = n (<^, -0?, h 

1 



+ (A;2 - t/ - sn(f/, fci) cn(f/, A;i) dn(f/, k{) ^ _ ^2g^2(^ fc^) 



and 



a^ 



2 RsiRa. — -Rp 



Rp{Ra. — R3 



. 2 ^ (/?a-i?3)(rn-fip) 

'''' ^ (i?^_i?p)(r„-i?3)' 
2 

(i?l-i?3)(i?2-i?4)(l-^2) 



{R± — R3)iRp — R±) 

E{ip, k) denotes the incomplete elliptic integral of the second kind, cn(-u, k) 
and dn(M, k) are the Jacobian elliptic functions, 

2 



J, 



2i?a 



(16) 
(17) 
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and Il{n,k) = n(7r/2, n, /c). J^, I±, Jr and Kr, in equations (|T2|), (|T4|) and 
(|15D are the integration constants. They depend on the number of turning 
points in r-coordinate between successive intersections with the equatorial 
plane (m) and on the sign of the radial velocity [sign(r)] at the intersections. 
For example, 

= K [mK{ki) - sign(r„+i) F{lp, k^)] . (18) 

Further details of the derivation are given in Appendix of Vokrouhlicky & 
Karas (1993) where the orbits under consideration are called the "Case II 
orbits." 



The above expressions can be considerably simplified in the special case of 
spherical orbits with r-coordinate constant. (It can be seen that r = const 
orbits do satisfy the geodesic equation due to integrability and separability 
of the geodesic motion in the Kerr spacetime; Chandrasekhar 1983). The 
gravitomagnetic precession of the spherical orbits around an extreme (a = 1) 
Kerr black hole was studied by Wilkins (1972). One can generalize his results 
to the case of spherical orbits r = Rs = const with an arbitrary value of the 
angular momentum parameter, | a |< 1. Taking into account equation (|^), 
we obtain 



{<i>n(-/i^,/i_^+) 



+a 



(£{r^ + a^) - a$) - S] K{fx_/fx+)} - 27r, (19) 



where A = — 2r + a^. The corresponding orbital period is 



Rs 



(r2 + a'^)S + 2ar {aS - $) 



(20) 



and the precession frequency is 



LT, 



r=Rs 




(21) 



Setting a = in equations (|T^)-(20) one arrives at 5</) = and P = 27iR^J^ 
which corresponds to Keplerian motion with vanishing nodal shift in the 
Schwarzschild metric, as expected. On the other hand, if a 7^ the dominant 
term in the asymptotic expansion of equation (O) for r — > 00 coincides 
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with the well-known result of Lense & Thirring (1918): Sep = 2TiaR~^^'^. The 
dependence of period P on the inclination fj, is only weak and equation (^) 
can be approximated by its value for equatorial orbits: 

P ^ 27r (r^/^ ± a) , (22) 

where +/— signs correspond to direct /retrograde orbits with respect to the 
rotating SBH. Analogously, the nodal shift of spherical polar orbits ($ = 0) 
can be obtained from equation (^) with /i_ = 1. Polar orbits were studied 
by Stoghianidis & Tsoubehs (1987). 

2.2. Motion of the pericenter 



In the present paper we are mainly interested in the frequency of the grav- 
itomagnetic precession. In the data, however, there will also be present a 
periodicity associated with another relativistic effect — the pericenter shift. 
We thus need to estimate the corresponding precession frequency Qp. In the 
Schwarzschild case 

Hp = 5(P/P, (23) 

where 

50 = -K{k) - 27r, (24) 

UJ 



IU2+Ui 



U3-U 



1 - 2ui 



i (1 - 2m) (1 - 2u2){-U (a', k') + -4^ 
1 Imi ^ ^ 1 — 2ui 

'■E{k') + (k'^ - a") K{k') 



+ [2a^k'^ + 2a 



2 a^- 



■ik'^)u(a',k') 



(25) 



Here, 



OJ = Vl — 2u2 — 4:Ui, 



2{u2 — Ui)uJ~ 



2 1 «2 

a =1 



k 



,2 U2 - Ui 



Ui Us - Ui 

where ui = l/Rs., U2 = 1/-Rp and = \ — ui — U2 are roots of the polynomial 
in M = 1 /r : 

U{u) = 2u^ -u^ + 2C-\ - (1 - S^)C-^ 
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[cf. "Case I orbits" in Appendix of Vokrouhlicky & Karas (1993)]. The 
angular momentum L is defined with respect to the axis perpendicular to 
the orbital plane of the star, not with respect to the axis of the disk like 
$ in the Kerr case. This definition is naturally more advantageous in the 
Schwarzschild case because the orbit is planar. 

An analogous effect of the pericenter shift can be expected in the Kerr case, 
too. Now, however, there is no unambiguous value of f2p because the orbit 
does not remain planar and the pericenter shift is mixed up with the gravito- 
magnetic precession. We calculate fip numerically by applying the mapping 
algorithm (j^) for the r-coordinate over a large number of revolutions of the 
star around the SBH. 

2.3. The mean value of the precession frequencies 

A specific value of the nodal shift 6(j) and corresponding precession frequency 
can only be associated with spherical orbits. In this case 50 is specified 
by equation (|19]) with three parameters — a, /i_ and Rg. The case of a quasi- 
elliptic orbit with nonzero eccentricity requires knowledge of four parameters. 
The appropriate parameters are a, pericenter distance Rp, and eccentric- 
ity e = (-Ra — -Rp)/(-Ra + -Rp)- The nodal shift oscillates between the maximum 
and the minimum values, 50max and 50min- Figures 1 and 2 show the graphs 
of 6(f){r) for some typical values of the parameters (here r denotes the radius 
of the intersection with the equatorial plane). For practical purposes one 
needs the mean value of the shift which characterizes an average taken over 
a large number of revolutions. We define 

(50) = 5(j){x) V{x) dx, (26) 

jRp 

where V{r) is the probability distribution for the radial coordinate of the in- 
tersections of the specified trajectory with the disk. V is normalized to unity, 
/ V{x) dx = 1. We apphed the following numerical procedure for evaluating 
equation (^): (i) at first, for a given set of orbital parameters {£, $, Q), we 
constructed the probability distribution by following the course of an ensem- 
ble of orbits with identical orbital parameters and randomly chosen initial 
conditions; (ii) then we applied equation (^ and performed the integration 
in equation (P^l). The resulting values of the mean nodal shift are tabulated 
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in the Appendix. One can see that the dispersion is often quite small and 
the mean value is a satisfactory approximation to the current value of 50, 
except for very close orbits. 

Analogously to (50), we introduce the mean orbital period (P) of an eccentric 
orbit and we calculate (f2p) by applying the mapping algorithm for the 
r-coordinate over a large number of revolutions of the star around the SBH. 

3. NUMERICAL MODELS 

3.1. Simple examples 

In this Section we will show how the relativistic precession becomes evident in 
a simulated signal from the source. We take into account general relativistic 
effects on photons coming from the source to the observer with no approx- 
imation. The method for integrating light trajectories in the Kerr metric 
has been discussed by many authors (Cunningham & Bardeen 1973; de Fe- 
lice, Nobili & Calvani 1974; Cunningham 1975 and 1976; Asaoka 1989). We 
employ an efficient approach described by Karas, Vokrouhlicky & Polnarev 
1992. We do not include the noise component which is of course present in 
a real signal. The reason is that the form of the noise depends on its partic- 
ular model and it does not affect the periodicities we are interested in. The 
way to account for the noise in the synthetic light curve is straightforward, 
in principle, if one adopts a specific model of its origin. Such a model can 
be based on mechanisms proposed by Moskalik & Sikora 1986, Chagelishvili, 
Lominadze & Rogava 1989, Abramowicz et al. 1991, or Baring 1992. 

The position of a distant observer is characterized by his inclination 6o with 
respect to the symmetry axis. We assume that the star-disk collisions mod- 
ulate the disk radiation at the moment when the star crashes through the 
disk from the far side to the near side with respect to the observer, i.e. once 
per revolution. The light travel time to the observer depends on the location 
of the intersection and it is thus periodically affected by the precession. Fre- 
quencies corresponding to the periodic modulation of the signal are evident 
in the Fourier transform of arrival times. 

The lensing effect enhances the radiation when the source is behind the black 
hole and thus contributes to the periodic modulation of the signal. As an 
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alternative to previously described Fourier transform of arrival times, one 
can detect relevant frequencies in the power spectrum of the photometric 
radiative flux. However, in this alternative approach the input signal for 
the Fourier transform depends on the model of the star-disk interaction (see 
below). We assume, for simplicity, that the shape of the observed signal 
from successive collisions has always an identical profile: a sharp onset and 
then an exponential decrease of the local luminosity in the static frame or, 
alternatively, in the disk co-rotating frame. In the latter case both the lensing 
and the Doppler effect contribute to the signal strength. Figure 3 is an 
example of typical light curves. We assumed a spherical orbit, r = 6, around 
a nearly-extreme Kerr black hole. It is not clear whether such a close orbit is 
astrophysically realistic because the process of capture is not well understood 
(Rees 1993). A more plausible configuration is treated in the next paragraph. 

The basic frequency in the power spectrum of the observed signal corresponds 
to the orbital motion. One can detect two fundamental lines in the power 
spectrum. The first line at higher frequency [Qo = Svr/P, see equation (|20|) ] 
corresponds to the orbital motion and the second one at lower frequency [^lt, 
see equation (pT])] corresponds to the gravitomagnetic precession. Figure 4 
shows the normalized power spectrum of the light curve from the previous 
Figure, strictly speaking the coefficient of spectral correlation as defined by 
Ferraz-Mello (1981). The power spectrum contains also linear combinations 
of fundamental frequencies. 

The shape of the light curve naturally depends on the particular model 
and position of the observer, however, we expect that identical periodicities 
caused by the precession of the orbit will be present in the signal indepen- 
dently of the details. The power spectrum constructed from arrival times 
has a simple form. In Figure 5 we applied the Fourier transform on time 
intervals elapsed between successive fiares in the light curve. Now, the shape 
of the spectrum is determined by orbital parameters of the star and angular 
momentum of the SBH but it depends only very weakly on other details of 
the model — decay time of the fiare, viewing angle of the observer, etc. This 
advantageous property is easy to understand because the arrival times of 
successive fiares are directly related to the precession of the orbit and we 
thus avoid the mo del- dependent features in the light curve. 
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As mentioned above, several effects tend to circularize the orbit of the star 
and for this reason spherical orbits are of special interest. However, the star 
can initially be captured in an eccentric orbit. Therefore, we also considered 
the case of quasi-elliptic orbits (e 7^ 0). The pericenter shift reveals itself as 
another line in the power spectrum. In Figure 6 one can detect lines corre- 
sponding to the orbital motion {Qo}- the gravitomagnctic precession {Qlt), 
and the pericenter shift (fip). The gravitomagnctic precession disappears in 
the Schwarzschild case and the power spectrum of an eccentric orbit is then 
rather similar to the case of a spherical orbit near a Kerr black hole. To dis- 
tinguish between the both extreme cases one needs an independent estimate 
of the radius and eccentricity of the orbit (Figure 7). 

3.2. The case of NGC 68I4 

At present there is no generally accepted explanation of the periodicity ob- 
served in the X-ray flux from the Seyfert galaxy NGC 6814. Among the most 
promising models are those which relate the periodicity to the orbital mo- 
tion of some object (unspecified for the time being) around the SBH. These 
models can in a natural way deal with the enormous stability of the period 
over several years. (For a discussion of different mechanisms which have been 
proposed in the literature see, e.g., Abramowicz 1992.) The model of a star 
colliding with an accretion disk has been considered as a viable model (Syer 
et al. 1991; Sikora Sz Begelman 1992; Rees 1993). The basic periodicity at 
approximately 12,200 s is interpreted as the orbital period of the star. It 
is a straightforward conclusion that if the SBH rotates then another peri- 
odicity associated with the gravitomagnctic effect should be present in the 
signal. This conclusion is independent of the unknown details of the star- 
disk interaction and it only assumes a nonzero viewing angle. Data which 
are currently available do not cover sufficiently long interval of time to reveal 
additional periodicities besides the basic one. For this reason we were un- 
able to give a definitive conclusion about the general relativistic precession 
in NGC 6814. This could be changed if data from ROSAT also show up the 
periodic component. However, already now one can see that the model of 
star-disk interactions faces serious difficulties (see below). Considering the 
time scale of the periodicity of NGC 6814 and the estimate of the black hole 
mass one can write 

P = 4.038 X 10-^M6(P), 
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f, 27r 
Plt = P- 



where P and Plt are measured in units of 12,200 s and Mq = M/ (lO^M©). 
Assuming e 0, P ^ 1 and Mq 1, one can estimate the radius of the orbit 
to be 

= (394 ^a)^/^ w 53.7. 

The periods associated with the gravitomagnetic precession and the peri- 
center shift are thus substantially longer than the orbital one, although one 
should remember that the estimate of the mass of the SBH in NGC 6814 and 
consequently the above value of Rg are rather uncertain (Bao 1992b); {6(f)) 
can be estimated from the tables which are given in Appendix. 

We constructed the predicted light curve in the way described in Fig. 3 and 
compared main frequencies in the power spectrum with those in observational 
data. The eccentricity was assumed to be limited to a very small value by 
the effects of orbital circularization; the low eccentricity is also required by 
the phase stability of the light curve, as discussed below. As a consequence, 
the dependence of Qp on e can be ignored. Also the inclination of the orbit 
is a free parameter but the dependence of all relevant frequencies on is 
very weak. One can see that the period related to the pericenter shift is in 
the range 16 P ^ Pp ;^ 20 P (the lower limit corresponds to a = 0, the upper 
limit corresponds to a — > 1). Time scales related to both the orbital motion 
and the pericenter shift are significantly shorter than the one related to the 
gravitomagnetic precession, Plt ^ 215P (indicated value is the lower limit 
corresponding to a — > 1). 

We considered the scenario in which the X-ray flares are associated with 
bright spots in the disk corona occurring in the place where the star crashes 
through the disk. After their creation the spots subsequently rotate with 
the disk matter, radiate with decreasing luminosity in their local frame, and 
eventually disappear after several revolutions. As the spots rotate, their ra- 
diation is periodically enhanced by the lensing effect which causes secondary 
maxima in the light curve. Secondary maxima can also be related to the event 
of the star crashing through the disk in the direction away from the observer. 
By comparing the EXOSAT (Mittaz & Branduardi-Raymont 1989; Fiore et 
al. 1992) and Ginga (Done et al. 1992) folded light curves, Abramowicz 
et al. (1993) demonstrated that the phase of important peaks in the light 
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curve remains remarkably stable. This fact imposes strong restrictions on 
the model and if it is confirmed one can conclude that the constant phase of 
individual patterns in the light curve can only be understood if the black hole 
does not rotate and the orbit of the star is nearly circular. The change in the 
shape of the light curve is then related to instabilities in the accretion process 
which affect the structure of the disk. A significant change in the accretion 
rate is evident from the change of the X-ray flux of NGC 6814 during the 
interval between the EXOSAT and Ginga observations. 

The peaks in the light curve are not regularly spaced. Within the framework 
of the discussed model it is very difficult to understand this non-regular (but 
constant in time) distribution of phases of the peaks in the light curve which 

has been reported at different ranges of energy. We remind that one of the 
peaks in the light curve is, according to the above model, associated with 
the creation of a bright spot at the moment of the star-disk collision. Sub- 
sequent peaks are gravitationally lensed images of the orbiting spot. These 
two mechanisms — star-disk interaction and lensing — are necessary because 
we deal with non-regular spacing of the peaks. The peaks originating by 
different mechanisms should very probably be distinguishable by comparing 
light curves at different energies. The distinction is not apparent in data, 
nevertheless, the model can definitively be ruled out only if the physics of 
X-ray flares in the disk corona is well understood; this subject is beyond the 
scope of the present paper. Further restrictions on the model come from the 
limits on P derived from the studies of the effects of gravitational radiation 
and star-disk collisions (King & Done 1993). Alternative explanations of the 
flares are based either on occultations of the central disk regions by matter 
ejected in the star-disk collisions or interaction of ejected matter with a jet. 
Within such models we need more data to restrict the angular momentum of 
the SBH and the eccentricity of the orbit, but these models face even more 
severe problems with the relatively narrow width of the profile of the peaks 
in the light curve and they require a precisely tuned inclination of the stellar 
orbit. 

4. CONCLUSIONS 

We considered the general relativistic precession of a star orbiting a super- 
massive black hole and colliding with an accretion disk in the nucleus of an 
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AGN. We derived formulae for the azimuthal shift due to the gravitomag- 
netic precession and perihehon shift during the free (geodesic) part of motion 
between subsequent interactions with the accretion disk. No restrictions on 
the orbital parameters and the black hole angular momentum were imposed 
except that they describe a stable bound trajectory around the Kerr black 
hole. Within the framework of this model our results restrict possible values 
of the angular momentum of the central black hole. To illustrate the pre- 
cession effect, we adopted a simplified model of the star-disk interaction and 
we determined relevant frequencies in the power spectrum of the observed 
signal. Both types of the relativistic precession which are relevant for our 
problem — the precession related to the pericentcr shift and the gravitomag- 
netic (Lense-Thirring) precession — expose themselves clearly in the power 
spectrum. We found the analysis of arrival times conceptually more trivial 
and at the same time more advantageous than the analysis of the complete 
photometric curve. 

Our conjecture is that the typical character of the power spectrum will remain 
conserved at some level in astrophysically more realistic models of the inter- 
action. Such an assumption is well-founded provided the orbital parameters 
of the star are not changed significantly during several periods associated 
with the precession motion {Qo ^ ^lt, ^p)- More realistic models must 
include the effects of tidal interactions, gravitational radiation and star-disk 
collisions on the orbital parameters (the work in preparation). 

We thank M. Abramowicz, C. Done and A. Lanza for very useful discussions 
concerning the problem of NGC 6814 and P. Haines for careful reading of 
the manuscript. We also thank the unknown referee and participants of the 
Relativity Seminar in Prague for helpful comments and suggestions which 
helped us to improve our work. 

APPENDIX 



We tabulate the azimuthal shift of orbital nodes due to the gravitomagnetic 
precession. The three tables correspond to a = 0.33 (Table 1), a = 0.67 
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(Table 2), and a = 1 (Table 3). In the tables, Rp denotes the pericenter 
distance in units of the gravitational radius e is the eccentricity of the 
orbit, and /i_ is the parameter characterizing inclination in the asymptotic 
region, r ^ _R+ [see equation (|^)]. The columns denoted by "+" and "— " 
correspond to direct ($ > 0) and retrograde ($ < 0) orbits, respectively. 
(The difference naturally disappears for polar orbits which are characterized 
by /i_ = 1 and $ = 0.) The nodal shift is given in the form 

IQ-^rad] , 

where (50) is the mean value of the shift and 50max5 ^4>min are the maximum 
and the minimum values of the shift (for details see the text.) In particular, 
for spherical orbits (r = const, e = 0) we obtain 50max = {^(p) = ^'Pmm, In this 
case, (Sep) is given by equation (|TP|). An ellipsis ". . . " in Table 3 excludes 
those combinations of parameters which do not correspond to a time-like 
geodesic. 
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FIGURE CAPTIONS 



FIG. 1. — The nodal shift 6(f){r) (in radians) per revolution for different eccen- 
trities of the orbit as a function of the radial distance of the intersection with 
the equatorial plane. Solid lines correspond to direct orbits, dashed lines 
correspond to retrograde orbits. Numbers given with each curve indicate the 
eccentricity. The inclination parameter is arccos /i_ = 45°, a = 0.9981. Two 
values of 6(j) for given r correspond to r > and r < 0, respectively. The 
case of spherical orbits is denoted by a star, " * " . 

FIG. 2. — As in Fig. 1 but for different inclinations and constant eccentricity. 
Now, e = 0.5 and the numbers indicate the value of arccos //_. In particular, 
arccos /i_ = 90° corresponds to equatorial orbits. Naturally, the curves of 
direct and retrograde polar orbits (arccos /i_ = 0) coincide. 

FIG. 3. — The numerical simulation of the observed light curve. Radiative flux 
is given in arbitrary units on the ordinate. The flux is periodically affected 
by the lensing effect (main peaks) when the source of radiation is behind 
the black hole. In this example we consider a source of radiation which is 
located in the equatorial plane in the place where the star intersected the disk. 
The position of the intersection is affected by the gravitomagnetic precession 
provided the central supermassive black hole rotates. As a consequence, 
the places of intersection are dragged in the azimuthal direction; the Figure 
covers 17 periods. The luminosity of the source is isotropic in the locally 
static frame and decreases exponentially with e-fold time r. The position of 
the observer at infinity is characterized by the viewing angle 6'o; 6'o = is the 
rotation axis of the black hole. The angular momentum parameter is chosen 
to be a = 0.9981 . The star moves in an inclined spherical orbit: r = 6, 
yU_ = 0.5 . The three cases correspond to {a) 6 = 80°, r = 0.2P, {b) 6 = 40°, 
r = 0.2P, {c) 6 = 80°, r = 1.5P where P\r=6 is the orbital period given by 
equation (^). Further details of the model are described in the text. 

FIG. 4. — Normalized power spectrum of the light curve from Fig. 3. The 
coefficient of spectral correlation is on ordinate and the spectral window on 
abscissa. 
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FIG. 5. — Normalized power spectrum of the arrival times for the signal from 
Fig. 3. The axes are as in Fig. 4. One can verify that the typical shape of 
the spectrum is quite insensitive to the details of configuration in a broad 
range of parameters. 

FIG. 6. — As in Fig. 5. The orbit of the star is now eccentric: e = 0.1, Rp = 6, 
H- = 0.5 . Parameters of the three cases (a)-(c) as before. 

FIG. 7. — Calculated light curves for the simple model described in § III a 

The case (a) is for the Kerr black hole with a = 0.9981 and a circular 
trajectory of the star at Rs = 40 (the value consistent with current estimates 
for NGC 6814). The long-term modulation of the light curve (three main 
peaks) is due to the gravitomagnetic precession (angular frequency Qur)- The 
case (6) is for a nonrotating black hole and an eccentric trajectory (e = 0.25, 
Rp = 40). Now the modulation is due to the pericenter shift and it has a 
different period which is given by (flp). 
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1. 3232 ±0 21«''_?-jl 2012.'_:;:^;' 1152 !_:;™| 



0. 1035 ±0 1246 ±0 794 ± 21 950 ± 51 638 it 29 760 ± 65 529 ± 31 628lg^ 450 di 31 532 it 66 
.25 1039 it 1243 it 797 it 21 948 it 51 640 it 29 758lgg 531 it 31 626+1^ 452 it 31 53llg5 

10 .5 1050 it 1233 it 806 it 22 941 it 49 647 it 30 752 i 62 536 it 33 62lt^5 456 it 32 527lg| 

.75 1074 it 1213 it 822 i 25 925 ± 45 659 it 33 740 it 57 546 ± 35 612^^^ 464 it 35 519l5g 

1. 1146 ±0 875 it 34 700 it 44 579 ± 46 491 ± 44 



0. 571 it 674 it 438 ±8 512 ± 16 351 ± 10 408 it 20 290 ± 11 335 ± 21 246 ± 11 283 it 20 
.25 573 it 673 it 439 ± 8 511 ± 15 352 ± 10 107 ± 20 291 ± 11 334 ± 20 246 ± 11 282 i 19 

15 .5 578 it 668 it 443 it 8 508 ± 15 355 it 11 404 i 19 293 it 12 332 it 20 248 ± 12 280 it 19 

.75 589 it 658 it 451 it 9 500 ± 14 361 it 12 398 it 18 298 it 12 328 it 19 252 it 12 276 it 18 

1. 624 it 476 it 11 380 ± 15 313 it 15 264 ± 15 



0. 375 it 435 it 287 it 4 330 it 7 230 ± 5 262 it 9 190 ±6 215 ± 9 160 ± 6 181 ± 9 
.25 376 it 434 it 288 it 4 330 it 7 230 ± 5 262 i 9 190 ±6 215 ± 9 161 ± 6 181 ± 9 

20 .5 379 it 431 it 290 it 4 328 it 7 232 ± 5 260 it 9 191 it 6 214 it 9 162 ± 6 180 it 9 

.75 385 it 425 it 295 it 4 323 it 6 235 it 6 257 it 8 194 it 6 211 it 9 164 ± 6 178 it 8 

1. 405 ±0 309 it 5 246 it 7 202 ± 7 171 ± 7 



0. 270 it 310 it 207 it 2 235 it 4 165 ± 3 186 it 5 136 ± 3 153 ± 5 115 ± 3 128 ± 5 
.25 271 it 309 it 207 it 2 235 it 4 166 ± 3 186 it 5 137 ± 3 153 ± 5 115 ± 3 128 ± 5 

25 .5 273 it 307 it 209 it 2 233 it 4 167 it 3 185 it 5 137 it 3 152 it 5 116 it 3 127 it 5 

.75 277 it 303 it 212 i 2 230 it 3 169 it 3 183 it 4 139 it 3 150 it 5 117 it 3 126 it 5 

1. 290 ±0 221 ±3 176 it 4 145 ± 4 122 ± 4 



0. 207 it 235 it 158 it 1 178 it 2 126 it 2 141 i 3 104 it 2 116 it 3 88 it 2 97 it 3 
.25 207 it 234 ±0 158 it 1 178 it 2 127 it 2 141 ± 3 104 it 2 115 it 3 88 it 2 97 it 3 

30 .5 208 it 233 it 159 it 1 177 it 2 127 it 2 140 it 3 105 it 2 115 it 3 88 it 2 96 it 3 

.75 211 it 230 it 162 it 1 175 it 2 129 it 2 139 it 3 106 it 2 114 it 3 89 it 2 95 it 3 

1. 221 it 168 it 2 134 it 3 110 it 3 92 it 2 



0. 165 it 186 it 126 it 1 141 i 1 101 it 1 112 i 2 83 it 1 91 it 2 70 it 1 77 it 2 
.25 165 it 185 it 126 it 1 140 it 1 101 it 1 111 i 2 83 it 1 91 it 2 70 it 1 77 it 2 

35 .5 166 it 184 it 127 it 1 140 it 1 101 it 1 111 it 2 83 ± 1 91 ± 2 70 it 1 76 it 2 

.75 168 it 182 it 128 it 1 138 it 1 102 i 1 110 it 2 84 i 1 90 it 2 71 it 1 75 it 2 

1. 175 it 133 it 1 106 it 2 87 it 2 73 it 2 



0. 135 it 151 it 103 it 1 115 it 1 83 it 1 91 it 1 68 it 1 74 it 1 57itl 62 i 1 
.25 136 it 151 it 104 it 1 115 it 1 83 it 1 91 it 1 68 it 1 74 it 1 57itl 62 i 1 

40 .5 136 ±0 150 it 104 it 1 114 it 1 83 it 1 90 it 1 68 ± 1 74 ± 1 58 ± 1 62 it 1 

.75 138 it 149 it 105 it 1 113 it 1 84 i 1 90 it 1 69 ± 1 73 ± 1 58 ± 1 62 i 1 

1. 143 it 109 it 1 87 it 1 71 it 1 60 it 1 



0. 114 it 127 it 87itl 96 ±1 69 it 1 76 it 1 57 ± 1 62 ± 1 48 ± 1 52 i 1 
.25 114 it 126 it 87itl 96 ± 1 69 it 1 76 it 1 57 ± 1 62 ± 1 48 ± 1 52 i 1 

45 .5 115 it 126 it 88 it 1 95 it 1 70 it 1 76 it 1 57 it 1 62 it 1 48 it 1 52 i 1 

.75 116 it 124 it 88 it 1 94 it 1 70 it 1 75 it 1 58 it 1 61 it 1 49 it 1 51 it 1 

1. 120 it 91 it 1 73 it 1 60 it 1 50 it 1 



0. 97 ± 108 it 74 it 82 ± 1 59 it 1 65 it 1 49 ± 1 53 di 1 41 di 1 44 i 1 
.25 98 ± 108 it 75 it 82 ± 1 59 it 1 65 it 1 49 ± 1 53 ± 1 41 ± 1 44 i 1 

50 .5 98 ± 107 it 75 it 81 ± 1 60 it 1 64 i 1 49 ± 1 53 ± 1 41 ± 1 44 i 1 

.75 99 ± 106 ± 76 ± 81 ± 1 60 ± 1 64 ± 1 50 ± 1 52 ± 1 42 ±1 14 ± 1 

1. 103 ± 7S ± (i2 ± 1 51 ± 1 l:-! ± 1 



TABLE 3 



Rp e a .2 A .6 .8 

i M- i + - + 2 ± Z ± Z ± Z. 

0. 8538 ±0 ... &548+l^ ... 53221^^® ... 4498+^^^ ... 3909^^^^ 
.25 8623 ±0 ... 6609^450 ... 5370tgi3 ■■■ 4537lgg| ... 39431^^3 

5 .5 8882 ±0 ... 6799lg3g ... 55201^^4 ... ^662+111 ... 4050l7g9 

.75 9341 ±0 ... 7151 ±723 ... 5808 ± 955 ... 4907^^°°^ ... 42651^^^ 

1. ... 8143i;™i e7S2tV,li 5752t??2* 5035 ± 2075 



0. 3066 ±0 4132 ±0 2366 ±81 3216 ± 480 1914tlio 2622 ± 586 16011^20 22041^^* 1373li2i 1897^646 
.25 3089 ±0 4121 ±0 2382 ± 84 3204^4^^ 19261^^^ 2611 ± 573 161ltj^^ 2194^^^^ 138lt}^J 1888^^^^ 

10 .5 3161 ± 4084 ± 2434 ± 93 3168 ± 434 19661^26 2576lv;33 16421"" 2163l;;33 1407l"g ISSgl^f,! 

.75 3300 ± 4005 ± 2535 ± 115 3094 ± 370 2044 ± 155 2508+4go ^™t{ll 2101^454 1460l^g2 18041*3^ 

1. 3695 ±0 2834 ±213 2284^273 1905^285 1630^276 



0. 1708 ±0 2276 ±0 1317 ± 31 1741 ± 130 1062 ± 43 1396i;g5 884+4g 1157lig" 754+4° 9831 
.25 1718 ±0 2269 ±0 1325 ± 32 1736 ± 128 1068 ± 44 1391^162 889^4^ 1153+ig4 758+|° 980+ 

15 .5 1752 ± 2245 ± 1349 ± 35 1716 ± 120 1086 ± 48 I375+J53 903 ± 52 1139^ "5 77015^ 968^ 

.75 1818 ±0 2195 ±0 1396 ± 41 1677 ± 105 1122 ± 56 13431i3b 932 ± 60 1112lj3g 793 ± 60 945^ 

1. 2018 ± 1542 ± 67 1235 ± 87 1023 ± 91 869lon 



0. 1129 ± 1476 ±0 870 ±16 1123 ±54 700 ± 22 895+gg 581 ± 24 738+71 494+23 623; 
.25 1135 ±0 1471 ±0 874 ± 16 1119 ± 53 703 ± 22 892^^^ 584 ± 25 7351^^5 496 ± 25 621; 

20 .5 1154 ±0 1455 ±0 888 ± 17 1107 ± 50 714 ± 24 883^^?, 592 ± 26 7271b7 503 ± 26 615; 

.75 1193 ±0 1423 ±0 915 ± 20 1083 ± 45 734 ± 27 864 ± 57 608 ± 29 712l^g 5I6I29 602; 

1. 1312 ± 1001 ± 30 800 ± 40 660 ± 42 559 ± 41 



0. 818 ±0 1051 ±0 630 ±9 798 ± 28 506 ± 13 634 ± 35 420 ± 14 521 ± 36 356 ± 14 439+3* 
.25 822 ± 1048 ±0 633 ± 9 796 ± 27 508 ± 13 632 ± 35 421 ± 14 520tp 358 ± 15 438^35 

25 .5 835 ±0 1037 ±0 642 ± 10 787 ± 26 515 ± 14 626 ± 33 427 ± 15 515 ± 34 362 ± 15 4341^4 

.75 860 ± 1014 ±0 660 ± 11 771 ± 23 529 ± 16 613 ± 30 437 ± 17 504 ± 31 370 ± 17 425I31 

1. 939 ±0 716 ±17 572 ±22 471 ± 23 398 ± 22 



0. 629 ±0 796 ±0 484 ± 6 603 ± 16 388 ± 8 479 ± 21 322 ±9 393 ± 21 273 ± 9 330 ± 20 
.25 632 ±0 793 ±0 486 ±6 602 ± 16 390 ± 9 477 ± 20 323 ± 9 392 ± 21 273 ± 9 329 ± 20 

30 .5 641 ±0 785 ±0 492 ± 6 596 ± 15 395 ± 9 473 ± 20 326 ± 10 388 ± 20 276 ± 10 327 ±19 

.75 659 ±0 769 ±0 505 ± 7 584 ± 14 404 ± 10 464 ± 18 334 ±11 381 ± 19 282 ± 11 321 ± 18 

1. 715 ±0 545 ± 10 434 ±14 357 ± 14 302 ± 15 



0. 503 ±0 628 ± 387 ± 4 476 ± 10 310 ±6 378 ± 13 257 ±6 309 ± 14 217 ±6 260 ± 13 
.25 505 ±0 627 ± 389 ± 4 475 ± 10 312 ±6 377 ± 13 258 ±6 309 ± 14 218 ± 7 259 ± 13 

35 .5 512 ±0 620 ± 393 ± 4 470 ± 10 315 ±6 373 ± 13 260 ±7 306 ± 13 220 ±7 257 ± 13 

.75 525 ±0 608 ± 403 ± 5 462 ±9 322 ± 7 366 ± 12 266 ±7 301 ± 12 225 ± 7 253 ± 12 

1. 567 ±0 432 ±7 344 ± 9 283 ± 10 239 ± 10 



0. 415 ±0 512 ±0 319 ±3 388 ± 7 256 ± 4 307 ± 9 211 ± 5 252 ± 9 179 ± 5 211 ±9 
.25 416 ±0 511 ±0 320 ±3 387 ± 7 256 ± 4 307 ± 9 212 ± 5 251 ± 9 179 ± 5 211 ±9 

40 .5 422 ±0 506 ± 324 ± 3 383 ±7 259 ± 4 304 ± 9 214 ±5 249 ± 9 181 ± 5 209 ± 9 

.75 432 ±0 496 ± 331 ± 4 377 ±6 264 ± 5 299 ± 8 218 ±5 245 ± 8 184 ± 5 206 ± 8 

1. 464 ±0 354 ± 5 282 ± 6 232 ± 7 195 ± 7 



0. 350 ±0 428 ±0 269 ± 2 324 ± 5 215 ± 3 256 ± 7 178 ±4 210 ± 7 150 ± 4 176 ± 7 
.25 351 ±0 426 ±0 270 ± 2 323 ± 5 216 ± 3 256 ± 7 178 ±4 209 ± 7 151 ± 4 176 ± 7 

45 .5 355 ±0 422 ±0 272 ± 2 320 ± 5 218 ± 3 254 ± 6 180 ± 4 208 ± 7 152 ± 4 174 ± 6 

.75 363 ±0 415 ±0 278 ± 3 315 ± 5 222 ± 4 250 ± 6 183 ± 4 204 ± 6 155 ± 4 172 ± 6 

1. 389 ± 297 ±4 236 ± 5 194 ± 5 163 ± 5 



0. 300 ±0 364 ±0 230 ±2 275 ± 4 185 ± 2 218 ± 5 152 ± 3 178 ± 5 129 ± 3 150 ± 5 
.25 301 ±0 363 ± 231 ± 2 275 ±4 185 ± 3 218 ± 5 153 ±3 178 ± 5 129 ± 3 149 ± 5 

50 .5 305 ±0 360 ±0 234 ± 2 272 ± 4 187 ± 3 216 ± 5 154 ± 3 177 ± 5 130 ± 3 148 ± 5 

.75 311 ±0 353 ±0 238 ± 2 268 ± 3 190 ± 3 212 ± 4 157 ± 3 174 ± 5 132 ± 3 146 ± 5 

1. 332 ±0 253 ±3 201 ± 4 165 ± 4 139 ±4 



